Effects of disorder on quantum fluctuations and superfluid density of a Bose-Einstein 

condensate in a two-dimensional optical lattice 
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We investigate a Bose-Einstein condensate trapped in a 2D optical lattice in the presence of weak 
disorder within the framework of the Bogoliubov theory. In particular, we analyze the combined 
effects of disorder and an optical lattice on quantum fluctuations and superfluid density of the 
BEC system. Accordingly, the analytical expressions of the ground state energy and quantum 
depletion of the system are obtained. Our results show that the lattice still induces a characteristic 
3D to ID crossover in the behavior of quantum fluctuations, despite the presence of weak disorder. 
Furthermore, we use the linear response theory to calculate the normal fluid density of the condensate 
induced by disorder. Our results in the 3D regime show that the combined presence of disorder 
and lattice induce a normal fluid density that asymptotically approaches 4/3 of the corresponding 
condensate depletion. Conditions for possible experimental realization of our scenario are also 
proposed. 
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I. INTRODUCTION 

Ultracold atoms in an optical lattice have opened a 
new theoretical and experimental window for investigat- 
ing fundamental problems related to condensed matter 
physics in a very versatile manner [H, [3, Q ■ Bose-Einstein 
condensates (BEC) trapped in an optical lattice allow for 
optimal control of the system's parameters: the depth of 
an optical lattice can be arbitrarily modified by chang- 
ing the intensities of laser beams Q, whereas the inter- 
atomic interactions can be controlled via the technology 
of Feshbach resonances [5(. Having at hand the possi- 
bility to shape the optical lattices and realize quasi-lD, 
quasi-2D and 3D BEC systems almost at will, therefore, 
an important direction of investigation consists in study- 
ing the properties of a BEC system along the dimensional 
crossover. 

Along this line, much work have been done [HQ. In 
particular, Orso et al.® theoretically showed that a 2D 
lattice can induce a characteristic 3D to ID crossover 
in the behavior of quantum fluctuations. Interest in 
this dimensional crossover needs to be renewed, how- 
ever, following the remarkable observation that disorder 
in quantum systems can have dramatic effects on the 
quantum fluctuations and superfluid density of a BEC 
0, 1, S, ES EH El 0, 0, El • For example, it has been 
observed that even a tiny amount of disorder in the con- 
fining fields leads to a fractioning of quasi- ID condensates 
in waveguide structures on atom chips [l6[ . Hence, a nat- 
ural question that immediately arises is how the disorder 
affects the 2D-lattice-induced dimensional crossover pre- 
dicted in Ref. 0. 

In this paper, we launch a systematic investigation on 



a BEC in the combined presence of 2D optical lattice and 
weak external randomness. The present work is divided 
into two parts. In the first part, we shed new light on the 
interplay between disorder and interatomic interaction 
along the 2D-lattice-induced dimensional crossover. Ac- 
cordingly, the analytical expressions of the ground state 
energy and quantum fluctuations are derived for a BEC 
in a tight 2D optical lattice in the presence of disorder 
using the Bogoliubov theory. Our results show that for 
a fixed atomic density, the BEC at small values of the 
lattice depth is in an anisotropic 3D regime. When the 
lattice depth increases, however, the BEC undergoes a 
dimensional crossover from an anisotropic 3D regime to 
a ID regime. The effects of disorder on such dimensional 
crossover is analyzed in detail. Our results for the ground 
state energy in the asymptotic ID regime without disor- 
der is in agreement with the exact Lieb-Linger solution 
expanded in the weak coupling regime We argue, 

therefore, that our result generalizes the exact ground 
state energy of Lieb-Linger model to that in the presence 
of weak disorder. All of our results in case of vanishing 
disorder are consistent with those obtained in Ref. 

In the second part of this work, we calculate the nor- 
mal fluid density of the condensate induced by disorder 
using the linear response theory. Accordingly, the trans- 
verse current-current response functions are calculated. 
These results in case of vanishing lattice depth recover 
corresponding ones in Ref. [181 ]. In particular, our re- 
sults in the 3D regime show that the combined presence 
of disorder and lattice induce a normal fluid density that 
asymptotically approaches 4/3 of the corresponding con- 
densate depletion [lj|]. 

The above scenario for disorder can be realized in cold 
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atomic systems using several methods in a controlled way. 
They include applying optical potentials created by laser 
speckles or multi-chromatic lattices 0, H3], introduc- 
ing impurity atoms in the sample [2l| and manipulating 
the collision between atoms J22|. The controllability of 
disorder in a BEC system makes the studies of disorder- 
induced effects extremely fascinating [23|, [U Hf| . 

The outline of this paper is as follows. In Sec. II, we 
introduce the Hamiltonian for a BEC trapped in a 2D 
optical lattice in the presence of weak disorder and then 
analyze this model using the Bogoliubov approximation. 
In Sec. Ill, the analytical expressions for the ground 
state energy and quantum depletion of such BEC system 
are obtained. Especially, we focus on the analysis of the 
combined effects of disorder and optical lattice on quan- 
tum fluctuations of this BEC system . Sec. IV presents 
calculations of reduction of the superfluid density due to 
disorder trapped in an optical lattice. Finally in Sec. 
V, we summarize our results and propose experimental 
conditions for realizing our scenario. 



II. HAMILTONIAN FOR A BEC IN BOTH 
PRESENCE OF A 2D OPTICAL LATTICE AND 
WEAK DISORDER 

The N-body Hamiltonian describing the Bose system 
has the form 



H — fj,N = / dr* f (r) 



ft 2 V 2 
2m 



- ti+Vopt(r) + V ran (r) 



+ !$t (r )$( r )j$( r)j (i) 

where ^(r) is the field operator for bosons with mass 
m, fi is the chemical potential, N = J dr^ (r)if! (r) is 
the number operator, and g — Airfi?a/m is the coupling 
constant with a being the 3D scattering length in the 
free space. In the Hamiltonian (JTJ), V ran (r) and V op t(v) 
respectively represent the external random potential and 
the 2D optical lattice. 

The 2D optical potential V op t(r) in the Hamiltonian 
([I]) is given by 



V opt (r) = sE R [sm 2 (q B x) + sin 2 (q B y)] 



(2) 



where s is a dimensionless factor labeled by the intensity 
of the laser beam and En — ti 2 q 2 3 /2m is the recoil energy 
with hqs being the Bragg momentum. The lattice period 
is fixed by n/d. Atoms are unconfined in the z-direction. 

Disorder in the Hamiltonian |T]) is produced by the ran- 
dom potential associated with quenched impurities (26j 



V ran (r) = 



E 



ri|), 



(3) 



with v(r) describing the two-body interaction between 
bosons and impurities, being the randomly distributed 
positions of impurities and Ni mp counting the number of 



r.;. To obtain the concrete form of the pair potential v(r), 
we need to investigate the scattering problem between a 
particle and a quenched impurity whose mass is taken to 
be infinite in the presence of 2D optical lattice. Here, we 
restrict ourself to the conditions of a dilute BEC system 
in the presence of a very small concentration of disorder. 
Thereby, the potential v(r) can be expressed by a pseudo- 
potential v(r) = gimp5(r)- Here, the gimp is the effective 
coupling constant that reads gimp = 2irh 2 b/m, where b 
represents the effective scattering length accounting for 
the presence of a 2D optical lattice [27} . 

In what follows, we will assume that the optical lat- 
tice is strong enough to create many separated wells that 
give rise to an array of condensates. Meanwhile, because 
of the quantum tunneling, the overlap between the wave 
functions of two consecutive wells are still sufficient to 
ensure full coherence even in the presence of disorder. 
In such case, one is allowed to use the Bogoliubov the- 
ory to study both equilibrium and dynamic behavior of 
the system at zero temperature. In addition, we also 
suppose that the chemical potential /i is small compared 
to the inter-band gap. By this assumption, we restrict 
ourselves to the lowest band, where the physics is gov- 
erned by the ratio between the chemical potential /i and 
the bandwidth 8t, where t is the tunneling rate between 
neighboring wells. Generally speaking, for fi <C St, the 
system retains an anisotropic 3D behavior, whereas for 
\x 3> 8t, the system undergoes a dimensional crossover to 
a ID regime. 

In the tight-binding approximation (28j . the low- 
est Bloch band of the BEC system can be writ- 
ten in terms of Wannier functions as 4>k x { x )<i ) k y {y), 



where <fik m ( x ) = J2i e w(x 
exp[-u 2 /2cr 2 ]/7r 1 / 4 C r 1 /2 wit h d/a 



Id) 



and w(x) = 
exp(-l/4 N /s). 

Expanding the field operators by the expression 'J'(x) = 
c i 'ke~ lkzZ 4 l k :r (x) 4>k {y)> the Hamiltonian {1} takes 
the form 



H = H - fiN 



2V 2- a k + qV-q a k' a k 
k,q,k' 



E^k'Kc-k'' 
k,k' 



(4) 



where e£ = h 2 k 2 /2m + 2t[2 — cos(k x d) — cos(k y d)] is the 
energy dispersion of the noninteracting model, V is the 
volume and g = Aith 2 a/m is an effective coupling con- 
stant with d = C 2 a and C — d /_f^ 2 w4 ( x )dx ~ dj \[2/ku. 
In Eq. , the Vk is the Fourier transform of V ran (r) 



Vk 



f 

V 



e lkT V ran (r)dv. 



(5) 



Here we choose that the randomness is uniformly dis- 
tributed with density n imp — N imp /V and Gaussian cor- 
related . Thereby, the two basic statistical properties 
of disorder are the average value (Vq) = gimpnimp and 
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the correlation function (V^V-k) — gf m p n imp/V . Here, 
the notation (...) stands for the ensemble average over all 
disorder configurations Q- 

In the following, we focus on the situation where the 
number of atoms in each tube is sufficiently large and 
refer to no as the condensate density. Under this as- 
sumption [2!|, we can neglect the Mott insulator phase 



transition which would occur only for extremely large 
values of the lattice depth. Throughout this paper, we 
consider the case when the strength of disorder is weak 
enough and Bogoliubov approximation still holds [26| . 
Hence, by proceeding in a standard fashion of Bogoli- 
ubov theory, the Hamiltonian (U) can be approximated 
as 



J 



Heff ~ " N 



V [ -fj.no + n V + -gn ) + 



k#0 



2 E ( £ k - M + 2 <?«o) (a k ak + a{ 



k^O 



k^O 



_ k O-k 



k^O 



(6) 



In the effective Hamiltonian ©, the only processes con- 
sidered are the annihilation of a pair {k, — k} into the 
condensate through the two-body interaction, the scat- 
tering of a single particle k into the condensate by disor- 
der, and the corresponding inverse processes. When the 
condensate becomes depleted, the last term in Hamilto- 
nian ([6]) is important [18, 30] and can be treated by mak- 
r t 1 2 

ing the replacement y Z)k/o "ik^k — * n 'Ek^o*k*ki 
where n' is a parameter to be determined later. We em- 
phasize here that the introduction of n' in Eq. ^ is 
equivalent to the introduction of an additional Lagrange 
multiplier as shown in the Refs. [1,0; EI- The effective 
Hamiltonian ^ is then diagonalized by the Bogoliubov 
transformation 



«k 



-"kCk - w k cl_ k 



Wkc k 



with 



WkC-k 



1 = 



/iVo(uk - VkY 

£k 



2 V El 



V k , 



V- 



(7) 



"03 _ -i 
E k 1 



E k = 



A) 2 



A), and A = g (no + n') — u. 



Consequently, the diagonalized Hamiltonian reads: 
H e ff - fiN = -Vun + Eg + 22 ^kC k c k , 

k^O 



(8) 



with the ground state energy E g reading 
1 



E„ 



^< 2 



nog 



Eu 



k^O 



N 



Tlimpgirnp 



n impgi mp 

V 



El 



(9) 



In conformity with general theorem [32j |. here we set 
A = g(no + n') — u — to ensure a gapless quasi- 
particle spectrum [3l|. This condition, together with 



n = no + l/^Sk^o^ik^k, determines n' and no as 
function of n. 



III. QUANTUM FLUCTUATIONS AND 
DIMENSIONAL CROSSOVER 



A. Beyond mean field correction to the ground 
state energy 

By replacing the sum in Eq. ^ with integrals in the 
continuum limit, we obtain 



v 



i 



1 gn 



2 Air a z \ff n o 



2/ 



Rgn 



AM^^ 9Q \^o 



2t 



,(10) 



with k = riimp/no and R = nb /a . The last two terms 



in Eq. (|10p are the beyond mean-field correction to the 
ground state energy due to disorder. The functions f(x) 
and Q(x) in Eq. (TIT))) are respectively defined as 



/(*) = 



ir r d 2 k 2^1(2, 2' 3 '-^(k)) 



2v/i J_ n (2tt)2 y^(k) 



(11) 



and 



Q(x) 



_ p (1 1 1 x~ \ 

7T / dk x 2^1V 2 ' 2' (l+x+xsin 2 (k x /2))' z > 



2- 



y^sin 2 (k x /2) + 1 



(12) 



In both Eqs. (jTTJ) and (fl2"|) , the variable x stands 
for 2t/gn that appears in Eq. (fTU|) . In Eq. (JTTJ , 
7(k) = 2 — cos(fc x ) — cos(ky). The function 2-F1 (a, b, c, rf) 
in Eqs. (fTTjl and (fT5|) is the hypergeometric function [33| 
and the integration over the transverse quasimomenta is 
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restricted to the first Brillouin zone, i.e. \k x \, \k y \ < 7T in 
Eq. (HJ) and \k x \ < n in Eq. (fT3j) . 

Eq. (|10p describes the ground state energy of a BEC 
in a combined presence of tight 2D periodic potential 
and weak disorder. Our model is characterized by three 
parameters: (i) the BEC parameter no, 3 , (ii) the concen- 
tration of disorder re = ni mp /no, and (iii) the ratio of 
effective scattering amplitudes b/a. The first one reflects 
the strength of interatomic interaction in the presence of 
periodic potential, whereas the other two are important 
parameters that describe the effect of disorder trapped 
in optical lattice. 



B. Lattice-induced dimensional crossover in the 
presence of weak disorder 

We now proceed to consider the asymptotic behavior 
of the ground state in the limit of x — > and x ^> 1, 
respectively. In the limit of x 1, corresponding to 
8t 3> fJ., the system retains an anisotropic 3D behavior; 
whereas for x — ► 0, corresponding to 8t <C n, the system 
undergoes a dimensional crossover to a ID regime. 

In the ID regime x — > 0, the f(x) saturates to the 
value 4-\/2/3 j^j. J n this limit, we can neglect the Bloch 
dispersion and Eq. d8j) approaches asymptotically to the 
ground state energy of a ID Bose gas in the presence of 
weak disorder 



L 



1 



9iDn 1D 



z9iDn 1D 




ink 



(n 1D g 1D f' 2 Q 



2f 



giDniD 



(13) 



with linear density ji\d = n$d 2 and coupling constant 
gio = g/d 2 . Here L is the length of the tube. In the case 
of vanishing disorder (re = 0), Eq. (|13p is in agreement 
with the exact Lieb-Linger solution of the ID model ex- 
panded in the weak coupling regime mgio/hniD -C 1 
QJ. We argue, therefore, that Eq. (TT3"|) with re ^ is 
the generalized Lieb-Linger solution of the ID model ex- 
panded in the weak coupling regime in the presence of 
weak disorder. 

In the opposite 3D regime x 3> 1, the functions 
(fTTj) and (TT2"|) approach respectively the asymptotic law 
f(x) ~ 1.43/ Vi- 16V2/(15nx) and Q(x) ~ -1/y/wx. 



Hence in the limit x 3> 1, Eq. (fTO)) takes the asymptotic 
form 
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(14) 



with m* = h 2 /2td 2 being the effective mass [28[ associ- 
ated with the band and a cr = —0.2Ady // m/m* being a 
further renormalization of the scattering length due to 
the optical lattice. [f|. The first two terms in Eq. (fT"4"]) 
are the mean-field contribution; whereas the remaining 
terms are generalized LHY correction in the presence of 
both optical lattice and disorder. In particular, the pres- 
ence of disorder gives rise to the last term in Eq. (| 14|) . 
Compared to the free case, the periodic potentials affect 
the ground state energy in two ways. First, the peri- 
odic potential increases the local atomic density, thereby 
enhancing interatomic interactions characterized by the 
effective scattering length a — C 2 a. Second, the lat- 
tice modifies the dispersion relation, which is captured 
by the effective mass to*. For vanishing disorder re = 0, 
Eq. (fl"4|) is reduced to the corresponding result in Ref. 
Q; whereas for vanishing optical lattice s = 0, our result 
recovers exactly the corresponding one in Ref. [26| . 



C. Quantum depletion 

The average number of particles with nonzero momen- 
tum represents a depletion of the condensate which can 



be calculated by N - N = Ek#o ( e k + 9 n o ~ E k) /2-Ek- 
In the continuum limit, the quantum depletion of a BEC 
trapped in a 2D optical lattice in the presence of disorder 
is calculated as 



N -No 



R 



N And 2 V no \gno ) ' Wirhd 2 \ no \gno ) ' 

where h(x) and K(x) are functions of the variable x = 
2t/gno that are respectively defined as 



h{x)-- 



d 2 k 



+ 00 



dt 



tt ( 2?r ) 2 JxjCk) y/t - X 7 (k) 



t + 1 

Vt 2 + 2t 



- 1 



and 



(16) 



K(x) = \ 



dk 



2-Pl(|,|,l,|i 



[l+x+xsin 2 (k/2)] 2 



2vr 



[l + x + xsin 2 (k/2)] ^Jx sin 2 (fc/2) 



(17) 
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The first term in Eq. (fT5|) arises from the two-body in- 
teractions of bosonic atoms trapped in the optical lattice. 
The second term, on the other hand, represents the effect 
of disorder on the quantum depletion. 

In the ID regime, h(x) ~ — ln(2.7a;)/v / 2 so the first 
term in Eq. (| 1 5[) diverges. It implies that in the absence 
of tunneling no real Bose-Einstein condensation exists 
which agrees with the general theorems in one dimension 
[30| . In the opposite 3D regime, the functions h(x) and 
K(x) respectively decay as 4/(37rv2a;) and 2/(nx). Ac- 
cordingly, the quantum depletion (115|) takes the asymp- 
totic form in the 3D regime: 



N-N, 







N 




(18) 



Eq. ([18]) generalizes the standard 3D result in the pres- 
ence of weak disorder in free space to that in the presence 
of an optical lattice. 



IV. REDUCTION OF SUPERFLUID DENSITY 
DUE TO COMBINED EFFECTS OF DISORDER 
AND OPTICAL LATTICE 

A. Longitudinal and transverse response functions 

In the second part of this paper, we apply the linear 
response theory to investigate the effects of disorder on 
the superfluid density of a BEC trapped in a 2D optical 
lattice. The general definition of the superfluid density is 
proposed by Hohenberg and Martin [34|. We emphasize 
that, unlike the condensate density, superfluidity is a ki- 
netic property of a system and the superfluid density is 
not an equilibrium quantity but a transport coefficient, 
and should be determined by the response of momentum 
density to an externally imposed velocity field [H, HE] • 
According to Kubo's formula [36| . the average momen- 
tum density (g(r, t)) induced by the external perturba- 
tion velocity field Sv(r ) is given by 



(ffi(r.*)> 



dt 



dr'Rij(rt,r t)Sv j (r)e et , (19) 



where 



9(r,t) 



2i 



* f (r)V*(r) - *(r)V$ f (r) 



(20) 



is the momentum density operator, and 



Rij(rt,rt ) = ( |3t(r,i),0j(r ,t ) 

dq duJ i[q.(iMr')-ai(tHf')l 

(2Vp 2^ 6 



%(q, W ),(21) 



is the momentum density correlation function, averaged 
with the ground state wavefunction for a BEC at rest 



without perturbation at zero temperature. The static 
current-current response function is given by [37| 



1 f + °° 



dui Rij (q, uj) 
2n uj — it 

(0\J^\n) (n\ Ji q |0 



(22) 



where w„ = (^n — £o)/^ and |0) is the ground state and 
the sum is performed over a complete set of excited states 
\n) with the energy ftw n . Here the current operator is 
defined by 



j x (y) 



2td 



^2 sin \ kx 



q X ( V ) 



a^k+q, (23) 



(24) 



In general, the superfluid density of a BEC distin- 
guishes between the low-frequency, long-wavelength lon- 
gitudinal and transverse responses of the system [18l. [30l . 

iSi, m, H3] 



X«(q) = ^^(q) + k-^)\/'q) 



The first term in Eq. (|25|) is the longitudinal component; 
it is parallel to q in both indices, i.e. J2i liiliQj 1 1 2 ) = 1j- 
The second term, the transverse component is perpendic- 
ular to q in both indices, i.e. qi{5ij — qiqj/q 2 ) = 0. 

A longitudinal probe corresponds to boosting the sys- 
tem and the entire fluid responds. On the other hand, 
a low-frequency transverse probe corresponds to a slow 
rotation of the system, which only couples to the nor- 
mal fluid while leaving the superfluid untouched. Hence 
the normal fluid density is defined in terms of the 
static transverse current-current correlation function as 

[H,[IJ HE [iffi 



lim xr(q). 

q^Q 



(26) 



The superfluid density is thereby defined as the difference 
between the total density p and the normal fluid density 



p s = p - lim xt (q)- 

q^O 



(27) 



There is an important point to stress here. By using the 
particle conservation, one can obtain p = limg_»oXi(q) 
[H, HE HE HE H3- This result is fixed by the model- 
independent /-sum rule [IE HE HE H3] • One is therefore 
easily led to directly write p s = lim 9 ^ [xl{<1) ~ Xr(q)]- 
This is not valid, however, for the calculation using the 
Bogoliubov approximation which violates the particle 
conservation by replacing a,Q with a c- number V No [l8| . 
In particular, computations involving the time evolution 
of zero-momentum particle states would be false. The 



6 



non-zero momentum states, on the other hand, are not af- 
fected. Thus for our present calculation of the superfluid 
density within the Bogoliubov' theoretical framework, it 
is the calculation of Eq. (|26p that can be trusted, not 
liniq^o xl (q) • Accordingly, the superfluid density can 
only be determined through Eq. (|27|) in the present con- 
text. 



B. Superfluid density 

The transverse response of a BEC is only due to the 
normal fluid, since the superfluid component can only 
participate in irrotational flow. For a BEC trapped in a 
2D optical lattice, the rotational symmetry is broken that 
gives rise to an anisotropic system. Consequently, the re- 
sponse in the unconfincd z-direction is different from that 
in the confined x-y plane. Nevertheless, the system is still 
isotropic in the x-y plane. Hence it's more convenient to 
consider a slow rotation with respect to the z-axis. We 
therefore take the transverse response function along the 
z-direction, given by 

^ h 2 k 2 f° 
Xzz (0) = 2^^||(|Vk| 2 > 
m Ez 

k K 

= ^pV^-oK^), (28) 

with I(x) being a function of a variable x = 2t/gno in 
the form 

r, \ r d 2 k 1 1 

T(x)= / - .„ - (29) 

> J-„ (2^ ^xW)T2 ^(k)T2) 2 

Two properties of the normal fluid density can immedi- 
ately be stated based on Eq. (|28|) : (i) the normal fluid 
component induced by disorder is not fixed by the long- 
wavelength properties of the excited quasi-particle. In 
order to assure the convergence of the integral in Eq. 
(128p , the behavior of the elementary excitation spectrum 
at high momenta is important; (ii) Eq. (|28p can be in- 
terpreted as the second-order term in the perturbation 
expansion of the normal fluid density in terms of weak 
disorder Vk- These two properties should be general and 
do not depend on the concrete form of disorder [l8j . 

The superfluid density is thereby found to be p s = 
p - p n with p n = Xzz{0). In the region x > 1, I(x) 
asymptotically approaches the function ^/2/(6ttx). Eq. 
(f2"8")l therefore takes the asymptotic form 

Xzz (0) = ^l^i?n (n a 3 ) 1/2 . (30) 
mi 

Eq. (f3"0")) generalizes the normal fluid density induced by 
disorder in the free space to that in the presence of opti- 
cal lattice. Accordingly, we conclude that in the asymp- 
totic anisotropic 3D regime, the disorder still generates 
an amount of normal fluid that is equal to 4/3 of the 



condensate depletion [18[ , despite the presence of optical 
lattice. 

To conclude this section, we note that in using Eq. 
(f2"T)) to determine the superfluid density, one has to spec- 
ify detailed probing approach and the limiting proce- 
dures thereof to distinguish transverse and longitudinal 
responses (l8l . [30L l37j . In rotationally invariant and uni- 
form systems, definition (|2"T|) has been a standard ex- 
pression to calculate the superfluid density [H, H3|. It 
would be less advantageous, however, for nonuniform sys- 
tems. There are other definitions for superfluid density 
0, |H, |3{|. Despite the formal difference, these var- 
ious definitions are fundamentally grounded in analyz- 
ing the linear response of a fluid to an external velocity 
boost. The underlying key quantity to determine su- 
perfluid components, therefore, is the momentum den- 
sity correlation function. In isotropic and uniform sys- 
tems, these definitions all lead to the same superfluid 
fraction within the linear response approximation. To 
determine the superfluid density in nonuniform systems, 
Refs. [1, [H, Hi| gave a more general expression. In the 
present context, both the definition (|27|) and the defi- 
nition in Refs. [H, l38l . [3Ql ] are appropriate. This is be- 
cause the presence of optical lattice has been accounted 
for through renormalizing the scattering length g and the 
kinetic energy term e£. [2g|. As a result, despite the pres- 
ence of lattice, one can study the quantum behaviors of 
a BEC as if the space is homogeneous [111 [29| . In partic- 
ular, in the limit 8t ^> p the effect of lattice is explicitly 
captured by the effective mass m* and the effective in- 
teraction g. Consequently in the asymptotic 3D region, 
a BEC trapped in an optical lattice can be effectively 
regarded as a uniform fluid composed of bosons with an 
effective mass m* and an effective coupling constant g 
[I, HI]. Thus within our present context, Eq. ^ still 
gives a good definition for the superfluid density partic- 
ularly in the asymptotic 3D region, despite the presence 
of lattice. 



V. POSSIBLE EXPERIMENTAL SCENARIOS 
AND CONCLUSION 

In this work, the physics of our model is captured by 
the interplay among three quantities: the strength of an 
optical lattice s, the interaction between bosonic atoms 
gn , and the strength of disorder R. All these quanti- 
ties are experimentally controllable using state-of-the-art 
technologies. The interatomic interaction can be con- 
trolled in a very versatile manner via the technology of 
Feshbach resonances. In the typical experiments to date, 
the values of ratio gno/En range from 0.02 to 1 [l], 0]. 
The depth of an optical lattice s can be changed from 
OEn to 32Er almost at will [1]. Disorder may be created 
in a repeatable way by introducing impurities in the sam- 
ple 12 1| . or using laser speckles and multi-chromatic lat- 
tices |19L[20. 23, 24]. Therefore, the phenomena discussed 
in this paper should be observable within the current ex- 
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perimental capability. 

We emphasize here that our investigation of a BEC 
in the presence of both optical lattice and weak disorder 
has been done within the Bogoliubov theory. Further 
improvement of the theoretical framework is to permit 
the treatment of the system properties for whole range 
of interatomic interaction strength, from zero to infinity, 
as well as arbitrarily strong disorder 0]. 

In summary, a BEC trapped in a 2D optical lattice 
in the presence of disorder is analytically investigated 
within the framework of the Bogoliubov theory. We 
focus on analyzing the combined effects of disorder 
and optical lattice on the quantum fluctuations and 
superfluid density of the BEC system. Accordingly, the 
analytical expressions of the ground state energy and 
the quantum depletion of the system are obtained. Our 
results show that the 2D lattice induces a characteristic 
3D to ID crossover in the behavior of quantum fluc- 
tuations. Furthermore, the effects of disorder on this 



dimensional crossover are discussed in detail. Especially, 
we argue that our results in the asymptotic ID regime 
is a generalized Lieb-Linger solution of the ID model 
expanded in the weak coupling regime in the presence of 
weak disorder. In addition, the normal fluid density of 
the condensate induced by the disorder is obtained by 
calculating the transverse response function. Our results 
show that in the 3D regime, the quantum depletion 
due to the combined effects of disorder and lattice will 
asymptotically approach 3/4 of the corresponding re- 
duction of the superfluid density. Finally, the conditions 
for possible experimental realization of our scenario are 
discussed. 
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